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We report on a significant discrepancy between recently published highly accurate variational 
calculations and precise measurements of the spectrum of Rydberg states in 8, Rb on the energy 
scale of fine splitting. Introducing a modified effective single-electron potential we determine the 
spectrum of the outermost bound electron from a standard WKB approach. Overall very good 
agreement with precise spectroscopic data is obtained. 
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I. INTRODUCTION 


reported in [5] choosing 


The spectrum of the outermost bound electron of an 
alkali atom like 87 Rb is hydrogen like, but lacks the n 2 - 
degeneracy of the eigenstates labeled by the principal 
quantum number n of the pure Coulomb potential mm 


E n ,i = — 


1 


(n — Si) 


2 ■ 


( 1 ) 


This effect is the well-known quantum defect Si, re¬ 
sulting from the interaction of the outermost electron 
with the ionic core of the atom and the nucleus. In 
a refined version of the statistical Thomas-Fermi the¬ 
ory 0 , an effective potential determining the interac¬ 
tion between the outermost electron and the nucleus can 
accurately be modeled by a spherically symmetric po¬ 
tential V eS (r; l) depending on the distance r from the 
center and depending on the orbital angular momentum 
l £ {0,1,2, ...,n— 1} 00,12: 


Ueff (r; l ) = -2 


Z e fi (r; l) 
r 


+ Vpoi (c 0 


( 2 ) 


Here the function Z e g (r; l) represents a position- 
dependent weight function that interpolates the value of 
the charge between unity for large r and charge number 
Z near to the nucleus for r —» 0, and V po i (r; /) represents 
a short-ranged interaction taking into account the static 
electric polarizability of the ionic core 0 El- 

Overall good agreement with spectroscopic data of al¬ 
kali atoms (but discarding the fine splitting) has been 


Z e fi (r;l) 
and 


1 + {Z-l)e~ ra 'W- re -™Ai) [ fl3 (/) + ra4 (z)] 

( 3 ) 


Vpoi (r;Z) 



( 4 ) 


A table of the parameters ai(Z), 02 (Z), 0 , 3 ( 1 ), 04 ( 1 ), a c , 
and r c ( l ) can be found in 0. 

In an attempt to also describe the fine splitting of the 
excitation spectrum of the outermost electron of 87 Rb, 
it has been suggested 0 to superimpose a posteriori a 
spin-orbit term 


Vso (' r]j,l ) 


^so (r\j,l) 
[l-a 2 U eff (r; l)] 2 ’ 


( 5 ) 


on the potential V e s (r; l), which then influences the spec¬ 
trum E n j i on the scale of fine splitting and the orbitals 
VVs,:/j( r ) accessible to the outermost electron. Here 


t r / ■ ;\ 2 1 (r, Z) .. . . 

Vso (r;j,l) = a - - - g(j,l), (6) 

r or 

and a = — ~ , denotes the fine-structure con- 

a b 137.036 

stant, and 


9 (j, Z) 


0 


if l = 0, 


2 


if l > 1, 


(7) 
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where j £ {Z — A /_|_ A}. To determine those orbitals 
(with principal quantum number n = n r + l + 1 and 
radial quantum number n r £ No), a normalizable solu¬ 
tion to the Schrodinger eigenvalue problem for the radial 
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wavefunction U n jj(r) = rR n jj (r) and associated eigen¬ 
values E n jj < 0 is required: 


d' 2 

dr 2 


where 


1(1 + 1 ) 

r 2 


+ V (r;j,l ) - E n j,, 


Un,j,l ( r ) = 0, 

( 8 ) 


V(r-,j,l) = V eS (r;l) + V S o(r-,j,l) (9) 


thus engenders a systematic (small) error of the matrix el¬ 
ements calculated in [7] on the fine-splitting scale. When 
employing substantially more STOs this error would cer¬ 
tainly become larger. With N = 500 STOs the discrep¬ 
ancy of these theoretical results with the high precision 
spectroscopic data, as shown in Tables |I] and [TTJ is far 
too large to be corrected by simply replacing Vso (r; j, l ) 
with Vgo (r; j, l). Hence another explanation is required. 


denotes the effective single-electron potential. 

A highly accurate variational calculation of the excita¬ 
tion spectrum of the outermost electron of 87 Rb has been 
carried out recently [7j, in which the authors expand the 
radial wavefunction of the Schrddinger eigenvalue prob¬ 
lem <[sj) in a basis spanned by 500 Slater-type orbitals 
(STOs). On the other hand, modern high precision spec¬ 
troscopy of Rydberg levels of 87 Rb has been conducted re¬ 
cently. Millimeter-wave spectroscopy employing selective 
field ionization allows for precise measurements of the 
energy differences between Rydberg levels [Sj. An inde¬ 
pendent approach is to perform purely optical measure¬ 
ments on absolute Rydberg level energies by observing 
electromagnetically induced transparency (EIT) [9, (TO]. 
However, there is a systematic discrepancy between vari¬ 
ational calculations and the spectroscopic measurements 
of the fine splitting 

A -E n ,l = E n l _i t - E n l+ i l (10) 

as shown in Tables U and El Given the fact that the er¬ 
ror bars of the independent experiments muni are below 
1.1 MHz down to 20 kHz, and on the other hand con¬ 
sidering the high accuracy of the numerical calculations 
presented in [7j, such a discrepancy between experiment 
and theory is indeed significant. 

So, what could be the reason for the reported discrep¬ 
ancies? First, it should be pointed out that in the vari¬ 
ational calculations [7] a slightly different potential was 
used, that is, 


V ( r;j,l ) = V eS (r; l) + Vso (r; j,l) ■ 


( 11 ) 


Certainly, within the first-order perturbation theory 
there exists no noticeable discrepancy in the spectrum of 
the outermost electron on the fine-splitting scale, when 
taking into account the spin-orbit forces with Vso (g j, 0 
instead of working with Vso This is due to the 

differences being negligible for r > Za 2 . However, since 
Vso (f; j , l) eventually dominates even the contribution of 
the centrifugal barrier term within the tiny region 

0 < r < a 2 Z , a subtle problem with a non-normalizable 
radial wavefunction U n j i(r) emerges when attempting to 
solve the Schrodinger eigenvalue problem for any l > 0 
with the potential Vso (r;j, l). Such a problem is absent 
when one works with Vso ( r \ j-, 0 0]- 


A variational calculation with the potential (11) em¬ 


ploying N = 500 normalizable STOs as basis functions 


II. QUASICLASSICAL APPROACH AND FINE 
SPLITTING OF THE HIGHLY EXCITED 87 RB 


In 1941 alkali atoms have already been studied in the 
context of modern quantum mechanics in the seminal 
work by Goeppert Mayer [3], who emphasized the excep¬ 
tional role of the l = 1 and l = 2 orbitals. According 
to Goeppert Mayer, the outermost electron of an alkali 
atom is governed by an effective r-dependent charge term 

Z eS (r) = l + (Z- 1 )F(r), (12) 

where the function F(r) has been determined by employ¬ 
ing the semi-classical statistical Thomas-Fermi approach 
to the many-electron-atom problem, posing the bound¬ 
ary conditions as lim r _>.o F(r) = 1 and lim^oo F (r) = 0. 
As discussed by Schwinger [T2], this approach ceases to 
be valid in the inner shell region Z~ x < r < Z~3 of 
the atom. Therefore, taking into account the fine split¬ 
ting in the spectrum of the outermost electron of alkali 
atoms a posteriori by simply adding the phenomenolog¬ 
ical spin-orbit term ^ to ([2]), resulting in the effective 
single-electron potential seems to be questionable on 
general grounds in that inner shell region. 

On a more fundamental level, the treatment of rela¬ 
tivistic effects in multi-electron-atom spectra requires an 
a priori microscopic description based on the well-known 
Breit-Pauli Hamiltonian mm 

Jl? = Jt? nI + J% s + Jt%. (13) 


Here is the ordinary nonrelativistic many-electron 
Hamiltonian, while the relativistic corrections are repre¬ 
sented by the perturbation operators and The 
perturbation term contains all the relativistic per¬ 
turbations like mass correction , one- and two-body Dar¬ 
win terms , and further the spin-spin contact and orbit- 
orbit terms, which all commute with the total angular 
momentum L and total spin S, thus effectuating only 
small shifts of the spectrum of the nonrelativistic Hamil¬ 
tonian Ji? nr . The perturbation operator Ji?f s on the other 
hand breaks the rotational symmetry. It consists of the 
standard nuclear spin-orbit , the spin-other-orbit, and the 
spin-spin dipole interaction terms, which all commute 
with J = L + S, but not with L or with S separately, thus 
inducing the fine splitting of the nonrelativistic spectrum. 

Although the proposed functional form of the poten¬ 
tial (11) is highly plausible on physical grounds out¬ 
side the inner core region r > Z~ s, prima facie it ap- 
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TABLE I. Fine splitting A E n ,i = \ for P states in [MHz]. 


State | n, l = 1} 

Exp. [11] 

Exp. [8] 

Theory [3 

Theory (this work) 

8P 

565.1(4) x 10 s 

N/A 

602.00 x 10 a 

567.75 x lO 13 

10P 

219.1(4) x IIP 

N/A 

231.87 x 10 a 

218.77 x lO 11 

30P 

N/A 

4246.30(5) 

4500.50 

4246.46 

35P 

N/A 

2566.41(32) 

2717.41 

2566.28 

45P 

N/A 

1144.09(13) 

1217.24 

1143.95 

55P 

N/A 

605.77(7) 

644.81 

605.68 

60P 

N/A 

460.76(5) 

480.32 

460.68 


TABLE II. Fine splitting AE n ,;=2 for D states in 

MHz], 

State | n, l = 2) 

Exp. [TT] 

Exp. gj 

Exp. QSI 

Theory [7[ 

Theory (this work) 

8D 

30.4(4) x 10 a 

N/A 

N/A 

113.17 x 10 a 

36.42 x 10 3 

10D 

14.9(2) x lO 11 

N/A 

N/A 

52.05 x 10 3 

16.56 x 10 a 

30D 

N/A 

452.42(18) 

452.5(11) 

1447.53 

456.13 

35D 

N/A 

279.65(10) 

280.4(11) 

894.84 

281.52 

45D 

N/A 

128.33(4) 

127.8(11) 

407.64 

128.98 

55D 

N/A 

69.17(2) 

69.4(11) 

223.71 

69.47 

57D 

N/A 

61.98(2) 

62.2(11) 

197.39 

62.24 


pears to be inconsistent to lump the aforementioned rela¬ 
tivistic many-body forces into an effective single-electron 
potential of the functional form (11), so that it pro¬ 
vides an accurate description also for small distances 
Z~ x <r<Z~\. 

In the absence of a better microscopic theory for an ef¬ 


fective single-electron potential V e s (r; j, Z) describing the 
fine splitting of the spectrum of the outermost electron 
in the alkali atoms, we introduce a cutoff at a distance 
r so (l) with Z~ x < r so (Z) < Z ~3 so that the effective 
single-electron potential is now described by the follow¬ 
ing modified potential: 


Vnod (V ji 0 


f Veff ( r ; Z) if 0 < r < r so (Z), 

[ Veff (r; l ) + V so (r; j , Z) if r > r so (Z). 


(14) 


The choice [2] 

r so (Z = 1) = 0.029483 x r c (Z = 1) = 0.0442825, 

(15) 

r so (Z = 2) = 0.051262 x r c (Z = 2) = 0.2495720, 


gives a surprisingly accurate description of the fine split¬ 
ting in the spectroscopic data for all principal quantum 
numbers n, see Tables [T] and [Til 

The calculation of the spectrum of the outermost 
bound electron is then reduced to solving the radial 
Schrodinger equation ([8]) with the modified potential 
Vnod (r; j, Z). The resulting spectrum is actually hydro¬ 
gen like, that is, 


E n ,j,l — 


( n - A jtl y 


(16) 


where Ay; denotes a quantum defect comprising also the 
fine splitting. In actual fact the quantum defect describes 
a reduction of the number of nodes n r of the radial wave- 
function for Z = 0,1, 2 as a result of the short-range in¬ 
teraction of the outermost electron with the ionic core 


of the atom. Because the higher the orbital angular mo¬ 
mentum quantum number Z, the lower the probability of 
the electron being located near to the center, it is clear 
that the quantum defect decreases rapidly with increas¬ 
ing orbital angular momentum Z. Therefore, Ay; is only 
notably different from zero for Z = 0,1, 2. 

Writing Ay; = Si + rjjj with r/y; <C Si, the fine splitting 
to leading order in a 2 is: 


A E, 


n,l 


= 2 


(n - 8{f 


(17) 


The quasiclassical momentum p = ff—Q of the bound 
electron with orbital angular momentum Z > 0, total an¬ 
gular momentum j = Z ± y, and taking into account the 

Langer shift Z(Z + 1) —» (Z + y) 2 in the centrifugal barrier 
Hang, is then given by 

(Z + 1 ) 2 ~ 

Q (r; j, Z, E) = V r2 2j + V mod (r; j, Z) - E. (18) 

For Z = 0 the centrifugal barrier term and the spin-orbit 
potential are absent. 
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FIG. 1. (Color online) The action integral v(j,l,E) associ¬ 
ated with the effective single-electron potential V mo d (r;j,l) 
vs. scaled energy for l = 0 (blue line), l = 1 (red line), 

l = 2 (green line), all for j = l + |. The curves for j = l — | 
only differ by a tiny shift proportional to a 2 . 



FIG. 2. (Color online) The quasiclassical momentum 
y/—Q ( r\j , l, E ) vs. scaled distance for l = 0 (dashed 
black), 1=1 (green), 1 = 2 (red), for E = E n jj correspond¬ 
ing to principal quantum number n = 57 and j = l + j . 
The main contribution to the quantum defect values in p4| 
originates from the inner core region r < r c (l). 


Considering high excitation energies E < 0 of the 
bound outermost electron, i.e. a principal quantum num¬ 
ber n> 1 , the respective positions of the turning points 
are given approximately by 


r (-) = 


( 1 + \) 


1 + V 1 + 0 + \) E 




-E 


1 + \ /1 + l 


E 


if l > 3, 


(19) 


if l > 1 , 


where 0 < l -C Of course for l = 0 only a single 

(large) turning point = Ye exists due to the absence 
of the centrifugal barrier. However, the lower turning 
points are strongly modified for l = 1,2 compared 
to the pure Coulomb potential case taking into account 
the core polarization. For l = 1,2 the relation A - ' 1 ( l ) ~ 
0.02 x r c ( l ) holds; that is, r* - ) (1 = 1) ~ 0.03472 and 
r(l = 2) ~ 0.12827 [2j. Since the cutoff r so (l) in (15) 
is substantially above those values of the lower turning 
points (/), a quasiclassical calculation of the fine-split 
spectrum of the bound outermost electron is reliable. 


For a chosen radial quantum number n r , the associated 
eigenvalues E = E rlJ j < 0 of the outermost electron now 
follow from the WKB patching condition HUES]: 


where v (j, l, E) denotes the action integral 


r ( + ) 


’ (j, l,E) = - f dry/-Q(r;j,l,E) 

d rp(r\j,l,E). (21) 


= tJ< 


Plotting the function v(j,l,E) versus 




for l = 


0 , 1,2 clearly reveals a linear dependence of the form 
v(j,l,E) = 7 =^ +c(j,/), see Fig. [I 
According to [ 6 |, for A, B , C, D SR, with A > 0, B > 
0, C > 0, and \D\ <C C the following equality holds: 




2B C D _ B /- BD 

+ V "" r* + ^ ~ 71 “ V + wVd 

( 22 ) 

For a pure Coulomb potential A = —E, B = 1, C = 
(l + \)~ and D = a 2 g(j,l). The corresponding action 
integral then reads 




v<&(3,1,E) = 


~ 0 + I) + 




*o+lr 


if l = 0 , 

if / > 0 . 
(23) 


It is thus found from WKB theory that the quan¬ 
tum defect associated with the single-electron potential 
kmod (?5 J; 0 is. 


! n r + 1 if l = 0 , 
n r + \ if l > 0 , 


( 20 ) 


A j i = lim 

E-> 0 - 


v(j,l,E)-v™ y >hE ) 


(24) 


Ignoring spin-orbit coupling, i.e. for a = 0 , one has 
Aj i = Si, the standard quantum defect. For l = 0 the 
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TABLE III. The values of quantum defect Aassociated with the Rydberg level n = 57 for l = 0,1,2. 


Quantum defect A 

Exp. [8] 

Exp. [101 

Theory [7] 

Theory (this work) 

Ai/2,0 

3.1312419(10) 

3.13125(2) 

3.12791 

3.13095 

Ai/2,1 

2.6549831(10) 

N/A 

2.65795 

2.65197 

A3/2U 

2.6417735(10) 

N/A 

2.64399 

2.63876 

Ai/2,1 — A3/2P 

0.0132096(14) 

N/A 

0.01396 

0.01321 

A 3/2,2 

1.3478971(4) 

1.34789(2) 

1.35145 

1.34851 

^ 5 / 2,2 

1.3462733(3) 

1.34626(2) 

1.34628 

1.34688 

A 3/2,2 — A 5 / 2,2 

0.0016238(5) 

0.00163(3) 

0.00517 

0.00163 


centrifugal barrier and the spin-orbit coupling term © 
are zero, so A jg —► Ai 0 = S 0 . 

The dependence of the quasiclassical momentum 
\J—Q (cX ^ E) on the scaled distance is shown for 
l = 0,1, 2 in Fig. [2j Clearly, it is the inner core re¬ 
gion (l) < r < r c (l) that provides the main con¬ 
tribution to the quantum defect values. We find, for 
l = 0, 2, that changing the fitting parameter a 3 ( l ) in 
([ 3 ]) from its tabulated value in [5] according to the 
scaling prescription a 3 (1=0) —t 0.814 x CI 3 (l = 0) and 
ci 3 (l = 2) —> 0.914 x 03 (l = 2), leads to a slight down¬ 
ward constant shift of the WKB-quantum defect. As a 
result of this change, the calculated WKB-quantum de¬ 
fect A; ± i ; then agrees well with the spectroscopic data, 
see Table III Such a change of a 3 ( l ) does not affect the 
fine splitting values A E n g though. We also find that the 
dependence of the fine splitting A E n g on the principal 
quantum number n is well described by for all n > 8, 
see Tables HI and UT1 

In actual fact, for r^\ which is a criterion that 

is always met for high excitation energies \J—E ~ 0 of the 
outermost electron, the uniform Langer-WKB wavefunc- 
tion (r) US:, 20, with A +s> considered as the only 

turning point, describes the numerical solution U n jg (r) 
to the radial differential equation © under the influ¬ 
ence of the effective modified single-electron potential 
(14) rather accurately }21| . Only very near to the sec¬ 
ond turning point at a distance smaller than r so ( l ), 
the Langer-WKB wavefunction (r) ceases to be 

a good approximation to the numerical solution U n jj (r) 


of the radial Schrodinger equation 


III. CONCLUSIONS 

In this work we reported a significant discrepancy be¬ 
tween experiment HUH! and highly accurate variational 
calculations [7J of the spectrum of Rydberg states of 87 Rb 
on the energy scale of the fine splitting. We discussed 
that the usual a posteriori adding of the relativistic spin- 
orbit potential to the effective single electron potential 
governing the outermost electron of alkali atoms is in¬ 
deed inconsistent inside the inner atomic core region. 
In the absence of a full microscopic theory that lumps 
all many-body interactions together with the relativistic 
corrections into an effective single-electron potential in 
a consistent manner, we suggested a modified effective 
single-electron potential, see (141, that enables a cor¬ 


rect description of the spectrum of Rydberg states on the 
fine splitting scale in terms of a simple WKB-action in¬ 
tegral for all principal quantum numbers n > 8. Modern 
precision spectroscopy of highly excited Rydberg states 
thus enables the probing of the multi-electron correlation 
problem of the ionic core of alkali atoms. This is cer¬ 
tainly a fascinating perspective for further experiments 
and theoretical studies. 
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